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1. Introduction 

Effective field theories (EFTs), which provide us a systematic perturbative scheme 
and a model-independent calculation method, have become a popular method to study 
hadronic reactions with and without external probes at low and intermediate energies. 
(See, e.g., Refs. [1, 2, 3, 4, 5] for reviews.) At very low energies, the Coulomb interaction 
becomes essential for the study of reactions involving charged particles. The first con- 
sideration of the Coulomb interaction in a pionless EFT was done by Kong and Ravndal 
(KR) for low energy S'-wave proton-proton (pp) scattering [6, 7]. They calculated the pp 
scattering amplitude up to next-to leading order (NLO). For loop calculations, they em- 
ployed dimensional regularization with minimum subtraction (MS) scheme and so called 
power divergence subtraction (PDS) scheme suggested by Kaplan, Savage and Wise [8, 9]. 
Then KR estimated a scattering length a(/i) for the pp scattering after separating off the 
Coulomb correction where [i is the scale for dimensional regularization. The leading or- 
der (LO) result of a(fi) was almost infinite at /j, — where m n is the pion mass [6]. 
In addition, the LO a(fi) was highly dependent on the value of fi. Including the NLO 
correction, they obtained a(/i = m^) = —29.9 fm [7] which is comparable to the value of 
the scattering length a np in the np channel, a np = —23.748 ± 0.009 fm 2 . 

The value of a(/x) deduced after separating the Coulomb and strong interactions is 
particularly important in the study of isospin breaking effects in S'-wave NN interac- 
tion [11, 12]. The accurate value of a np is well known as quoted above, while the values 
of the scattering length in the nn channel (a nn ) and in the pp channel (a pp ) still have 
considerable uncertainties. 

There exists no direct nn scattering experiment because of the lack of free neutron 
target. The values of a nn have been deduced from the experimental data of 7r~d — > nrry 
and nd — > nnp reactions. Recent publications suggest a nn = —18.50 ± 0.05(stat.) ± 
0A4(syst.) ± 0.30(t/i.) fm from the n~d — > nn^ process [13] and a nn = —18.7 ± 0.6 
fm [14], —16.06 ± 0.35 fm [15] and —16.5 ± 0.9 fm [16] from the nd — > nnp process. As 
seen, the values of a nn have significant errors compared to that of a np , and the center 
values do not seem to converge yet. 3 

For the pp channel, a very accurate value of the scattering length ac = —7.828 ± 0.008 
fm [19] and ac = —7.8149 ±0.0029 fm [20] are available from the low energy pp scattering 
data. It contains however contributions from both strong and electromagnetic interac- 
tions, and thus we need to disentangle the strong interaction from the electromagnetic 
interaction. It was pointed out in potential model calculations that there is a consid- 
erable model dependence in deducing the value of the strong scattering length a pp from 
ac [19, 21]. Some literature shows a pp = —17.1 ±0.2 fm [19], while a heavy-baryon chiral 
perturbation theory results in a pp = —17.51 ~ —16.96 fm [22] with uncertainties slightly 
larger than those from the potential models. 

In this work, we employ the pionless EFT [23] including the Coulomb interaction 
between two protons [6, 7] and calculate the pp scattering amplitude with the strong 

2 See, e.g., Table VIII in Ref. [10]. 

3 Recently, there were proposals to determine the value of a nn more precisely by employing a formalism 
of EFT, from the n~d — > nnj reaction[17] and neutron-neutron fusion, nn — > de~D e {18]. 
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NN interactions up to next-to-next-to leading order (NNLO). Our main motivation of 
this study is to see how the value of strong scattering length a(/i = m^) = —29.9 fm 
obtained by KR from NLO calculations may be improved by the inclusion of a higher 
order correction. We find that the NNLO corrections turn out to be quite small but there 
is a considerable dependence of the scattering length a(/j,) on the renormalization schemes 
and the scale parameter \i. 

This paper is organized as follows. In Sec. 2 we briefly review the effective range 
formalism for the pp scattering. In Sec. 3 the pionless strong effective Lagrangian up to 
NNLO is introduced. In Sec. 4 we calculate the S'-wave pp scattering amplitude up to 
NNLO. In Sec. 5, we discuss regularization method and renormalization schemes employed 
in this work. We renormalize low energy constants (LECs) that appear in the strong NN 
interaction up to NNLO by effective range parameters employing MS-bar (MS) and PDS 
schemes and obtain numerical results for the strong scattering length a(/x) and strong 
effective range r(/x). Discussion and conclusions are given in Sec. 6. In Appendix A we 
show detailed expressions of the amplitudes in NNLO. Detailed calculations of the loop 
functions employing the dimensional regularization and MS and PDS schemes are given 
in Appendix B. 

2. Proton-proton scattering in effective range theory 

The amplitude of the pp scattering can be decomposed as [24] 

T = T C + T SC , (1) 

where T c is the pure Coulomb part and T S c is the "modified" strong amplitude whose 
S'-wave channel we calculate up to NNLO in pionless EFT below. 

The incoming and outgoing scattering states | x 3/^ ± ' > ) with the potential V — Vc + Vs 
where Vc and Vs are the Coulomb and strong potentials, respectively, are represented in 
terms of the Coulomb states as 

oo 

\*? ) ) = T,(G£ ) Vs) n \1>? ) ), (2) 



?1=0 



where G^ is the incoming and outgoing Green's function 

Gg° (E) = . (3) 

E — H — Vc±ie 

Here Ho — p 2 /M is the free Hamiltonian of two protons and Vc = e 2 / (Aixr) is the repulsive 
Coulomb potential. The Coulomb state 1^4^) is obtained by solving the Schrodinger 
equation (H — E)\ip^} = with H = H + V c and thus one has 

= V + d c ] Vc] \P) , (4) 
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where \p) is the free wave state. The normalization of is such that (V^IV^) 

(27r) 3 S ( - 3 \p — q). The amplitude T S c is thus obtained by 



Tsc(i?,P) = E(4 _) l^(^ +) V5)l4 +j )- 



-,(+)\ 



(5) 



n=0 



For / = state one has the amplitude 



T, 



1=0 

sc 



47T 



M p cot 5q — ip ' 



(6) 



where 07 is the Coulomb phase shift 07 = argT(l + I + irj) with rj = aM/(2p). In the 
effective range expansion with the Coulomb interaction, the modified strong phase shift 
5i for / = in low energy pp scattering is represented by effective range parameters [25] : 



Cl p cot 5 + aMh(r)) = — - + \r p 2 - Prgp 4 + 

dc 2 



where C 2 T] = 2vr??/(e 2 ™' - 1) and 



h(i]) = Re ip(irj) — In r\ . 



(7) 



(8) 



■^-function is the logarithmic derivative of the Gamma function and Re ip(irf) = rf Yf^Li v^+rj 2 ) ' 
Ce] Ce is the Euler's constant, Ce = 0.577215 ■ • •. Effective range parameters ac, r , P 
are modified scattering length, effective range, effective volume, respectively. 

3. Effective Lagrangian 

Pionless effective Lagrangian for strong 5-wave NN interaction up to NNLO reads [23, 

26] 



2m 



N-C 



N 



4° 2 



<-+2 it 

N^pCso) D N N T P es ) N + h c 



(9) 



where is the covariant derivative, D= \{D — D), and P^ So ^ is a projection operator 
for the two-nucleon 1 S states, P^ So ^ = ^gO^Ta- Note that we retain two low energy 

constants, C 4 and (7 4 , in NNLO. 

The strong NN potential is expanded in terms of small momentum as 



V s = V + V 2 + V A + ■ ■ ■ , 



(10) 
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where V , V 2 , Va are LO, NLO, NNLO potential, respectively, and the matrix elements of 
them are obtained from the Lagrangian in Eq. (9) as 

(q]V \k) = Co, (11) 
(q]V 2 \k) = }-C 2 (f + k 2 ), (12) 

(q]V 4 \k) = ^fk' + ^C^ + k 4 ), (13) 

— * 

where \q) and \k) are the intermediate free two-nucleon outgoing and incoming states, 

— * 

respectively: 2q and 2k are the relative momenta for the two protons. 

In this work we employ the standard counting rules of the strong NN interaction with 
the PDS scheme in Refs. [7, 8]. (We will discuss the PDS scheme in detail later.) For 
the strong potential, the LO term Co is counted as Q^ 1 order, where Q denotes the small 
expansion parameter, and is summed up to an infinite order. The NLO (C2) and NNLO 
(C4, C4) terms are counted as Q 2 and Q 4 , respectively, and expanded perturbatively. 4 
We treat the Coulomb interaction non-perturbatively using the Green's function in 
Eq. (3). We do not include higher order QED corrections such as the vacuum polarization 
effects reported in Refs. [27]. 

4. Amplitudes 

The amplitude T^f for the «S-wave pp scattering can be written as 

T l s f = T^ + T^ + T^ + .-, (14) 

where T^, T<£l, T$l are LO, NLO, NNLO amplitudes, respectively. By inserting the 
strong LO potential V in Eq. (11) into the amplitude T sc in Eq. (5), we obtain the LO 
amplitude Tg^ in terms of loop functions ip and J : 

42 = t^iWcmri^) = i c _f^ {p) , (15) 



n=0 



*•<» = f^to'fw**®' (16) 



where 



Detailed calculations for the functions ipo and Jo are given in Appendix B. T§q is sum- 
mation of the LO strong potential Vo, that is, the C terms summed up to the infinite 
order. 



4 Note that by changing the LECs C4 and C4 in another linear combination, e.g., C4 = C 4 + C' 4 and 
Ca = C 4 — C4, one can easily see that the term proportional to C' 4 in Eq. (13) vanishes when \q\ = \k\. 
The C 4 term, so called off-shell term, is redundant and vanishes when the external legs of the potential 
go on mass-shell. 
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(a) 




(b) 





(c) 





(d) 






Figure 1: NLO diagrams for the S"-wave pp scattering. Gray blobs denote the two- 
proton Coulomb Green's function Gq , and two nucleon contact vertices denote the strong 
potential: the (black) circle and the (red) square represent LO (Co) and NLO (C2) vertices, 
respectively. Small double dots stand for the summation of Co terms up to the infinite 
order. 



At NLO we have four diagrams shown in Fig. I. 5 They are proportional to C2 coming 
from V2, whereas the Co terms are summed up to the infinite order. The NLO amplitude 
is written in terms of the loop functions ipo, ip2, Jo and J2 as 



with 



ip 2 (p) 

Up) 



r (2,a-d) 
SC 



d 3 k 



C 2 ^ 



;i - c j ) 



[lp2 + C (^ J 2 ~ 1P2J0)} 



d 3 k 



(2 



7T 



d 3 q d s c? j 2 _ ~ (+) 



(2tt) 3 (2 



7T 



d 3 q d 3 <? 
(2tt) 3 (2tt) 3 



->2 



(18) 

(19) 
(20) 



(2) 

Details for ip2 and J2 are given in Appendix B. The NLO amplitude consists of one C2 
and a summation of the Co terms up to the infinite order. These LO and NLO amplitudes 
have already been obtained by KR in Ref. [7]. 

At NNLO we have three sets of diagrams shown in Figs. 2, 3, and 4. From the first and 



'Figures were prepared using the program JaxoDraw [28] provided by L. Theussl. 
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Figure 2: Set 1 of NNLO diagrams. See the caption of Fig. 1 for details. 



second sets of diagrams shown in Figs. 2 and 3, respectively, we see two NLO corrections 
to the amplitude and thus the NNLO amplitudes obtained from the first and second sets of 
diagrams in Figs. 2 and 3 are proportional to Cf. The NNLO amplitudes corresponding 
to the diagrams in Fig. 2 can be written in terms of the functions ipo, ip%, Jo and 
whereas to express the amplitudes for the diagrams in Fig. 3 we need a new function J 2 2 
given below. In the third set of diagrams shown in Fig. 4, we have one NNLO correction 
to the amplitude and the NNLO amplitudes for the diagrams in Fig. 4 are proportional 
to C4 or C4. Explicit expressions of the NNLO amplitude from each of the diagrams are 
given in terms of ipi with i — 0, 2, 4 and Jj with j = 0, 2, 22, 4 in Appendix A. 

Summing up the amplitudes obtained from the diagrams (a) to (h) in Figs. 2 and 3 we 
have 

Tfr h) = 4(1 Jo) 3 {^(1 - Q) Jo) + i>\ Ml - CoJof 

+2VWW1 -ClJl) +iJ 2 J2(C J 2 )(3 + CoJo)} , (21) 

where 

*■ = lSwr^ wm ^ f ' (22) 

whose details are given in Appendix B. 
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Figure 3: Set 2 of NNLO diagrams. See the caption of Fig. 1 for details. 



Summing up the amplitudes for the diagrams (i) to (1) in Fig. 4 gives us 



T. 



(4,i-0 



1 



Ca 



so 



2 (1 - Co Jo) 2 
1 C 4 



2 (1 - C J ) 2 



^(i - c j ) 2 + 2^ 2 c j 2 (i - c j ) + r ci Ji 



[-04 + C O (-0 O J4 - -04 Jo)] IpO 



where 



V>4 
J 4 



d 3 k ,(->, ^ r4 _ /" ^ T4 /(+)/ 



d 3 <f d 3 ^ ^ \M+)\-\ f rf3( ? <Pf 

-,q w \Gc w = 



(2vr) 3 (27r) 3 ^ ^ 10 J (2vr) 3 (2tt) 3 

Calculations of ^4 and J4 are given in Appendix B. 



(23) 

(24) 
(25) 



5. Regularization method and renormalization schemes 

In the calculation of the loop functions J , J2, J22 and J 4 in Eqs. (17), (20), (22), (25), 
we encounter infinities and employ the dimensional regularization. We also employ the 
PDS scheme, suggested by Kaplan, Savage and Wise [8, 9], in which one subtracts the 
poles in d = 3 as well as those in d — 4 space-time dimensions so that one obtains an 
expected perturbation series in the expansion of the NN potential in Eq. (10) with a 
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Figure 4: Set 3 of NNLO diagrams. Two-proton contact vertices represented by (blue) 
diamonds denote strong NNLO potential V4. See the caption of Fig. 1 for details. 



given scale /1 of the theory. We may check the convergence radius, e.g., for the C2 term 
(relative to the C term) in Eq. (10) and have A 20 (/i) = ^C (^)/C 2 (n) = 147 (30.6) MeV 
with (without) the PDS terms at p = m n . Thus a formal convergence of the perturbative 
series of the NN potential in Eq. (10) is improved thanks to the PDS term, and the 
theory would be valid up to p ~ A 2 o — 140 MeV, which is the large scale we assumed in 
the pionless theory. 

The loop functions can be decomposed into a finite term and an infinite one, e.g. 
Jo = Jo™ + Jq w with Jg m = —^—Hirj) (the definition of the H{rj) function is given in 
Appendix B) and 



J; 



div 



M 
An 



aM 2 



li 



57T 



3C E + 2 + In 



7r/r 
a 2 M 2 



(26) 



where Jq %v 



is calculated by the dimensional regularization in d = 4 — 2e dimensions and 
the PDS scheme. The first term proportional to the scale \x in the r.h.s. of Eq. (26) is the 
PDS term and Ce is the Euler's constant mentioned earlier. The scattering amplitudes 
should be identical after renormalization even if another renormalization scheme such as 
off-shell momentum subtraction scheme discussed in Refs. [26, 29] is employed. However, 
a(fjj) and r(/x) do depend on the renormalization schemes along with the value of the 
renormalization scale /x. So, to be consistent with KR, we calculate all the loop functions 
Jj with i — 0, 2, 22, 4 and the wavefunctions ipj with j = 0, 2, 4 by using the dimensional 
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regularization and the PDS scheme in Appendix B. 

The .S-wave pp scattering amplitude in terms of the effective range parameters is given 

by 



rpl = 

1 SC 



4.7T 



fi2 2ia 



M -aMH(rj) - ± + \r Q p 2 - Prjjp 4 + 



(27) 



and thus one has 



— — + \r Q p 2 - Prlp 4 + 
ac 2 



= aMH(rj) 



47T ijj, 



M T f) 



Ait 1 



P (2) 

■sc 



T. 



(4) 



SC 



(0) 



sc 



t: 



(0) 



+ 



■sc 



sc 



T, 



(0) 



sc, 



sc 



+ 



(28) 



Comparing the coefficients of the terms proportional to p°, p 2 and p 4 in both sides of 
Eq. (28), we have 



1 

a c 



1 

+ 2 r ° 



-Pri = 



M\C J ° + C% 

1 Ca Co 



aMfj, + \{uM) 2 + C ^-(aM) 2 /i 
2 48 



+ o(« 3 



9i 



-2 



2Cl + 3Cl 




aM)n 



+ 0{a 2 ) , 



IC4 IC4 
,2C| + 2Cf 



(29) 
(30) 
(31) 



where we have expanded the r.h.s. of Eqs. (29) and (30) in the order of the fine structure 
constant a and neglected the a 3 (a 2 ) and higher order terms in Eq. (29) (Eq.(30)). With 
three effective range parameters, we cannot determine the four LECs uniquely. There are 
some arguments which can constrain the values of Ca and Ca. The C4 contribution in 
Eq. (29) is of the order of /i 2 , and thus the first C4 contribution term is of the lower order 
of fi than the C4 term. 6 For this reason, the C4 term is treated as an order higher than 
the C4 one [23], and consequently the C4 term does not appear (at NNLO) in Eq. (29). 
The other argument is based on the offshell-ness of a term proportional to C4 — C4 [26]. 7 
In this case, the term proportional to C 4 — C4 is redundant and thus can be removed 
by assuming C4 = C±. Because both arguments seem to have some grounds, to check 
the dependency of the results on the values of C4 and C4 we consider the three cases: 1) 
Ca = (Ref. [23]), 2) C 4 = Ca (Ref. [26]), and 3) C 4 = 0. 

In Eq. (29) there is the Jq iv term explicitly given in Eq. (26). In the MS scheme used 



by KR [6, 7] one subtracts the infinite term 



M 2 1 



8tt 



from the J$ w . 



One can use another 



6 Note that fi is regarded as a large scale, i.e., /i — 
7 See the footnote 4. 
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scheme called MS scheme, in which finite terms are subtracted together with the infinite 



term so that 



aM 2 
8tt 



\ - C E + 1h(4tt) 
M 



J, 



MS 



All 



is subtracted. Then we have 
aM 2 



An 



In 



// 



2aM 



+ 1-C 



E 



(32) 



This leads to a significant subtraction scheme dependence in the scattering length a(/i). 

6. Numerical results 

We may define the strong scattering length and the effective range, respectively, in the 
zeroth order of a as [7] 



1 



47T 



1 AnC 2 (fx) 



(33) 



a(n) M Co(pl) 

Inserting the expressions of a(/i) and r(fj) in Eqs. (33) into Eqs. (29) and (30), we have 



J LO 



a{n) 
r - (aM) 



+ 



a(fi) 



+ 



NLO 



D 3 Pr*n + L> 4 - 



a(fx) 



J NNLO 



where 



a(/x) 
1 



LO 



a c 



+ aM 



In 



2aM 



ac ^ 



+ 1-C 



(34) 
(35) 



E 



a(/x) 
1 



NLO 



--aMr fi - (aMf 



r /i 



1 7T : 

4 ro + r 2 ^ 

ac 



a(/i) 



(aMf 



NNLO 



12 —-li 

ac ^ 



(36) 
(37) 
(38) 



and the term linear in aM in Eq. (35) is the NNLO correction to r(/i). We have three 

{Di, D 2 , D 3 , -D4}, because of the additional constraints 



set of coefficients, -^ x (=i,2,3) 
imposed on the LECs C4 and C4 mentioned before Eq. (32). X\ = {1, 0, 4, 0} corresponds 
to the case 1) C4 = 0, X 2 = {7/6, 1, 10/3, 1/6} corresponds to the case 2) C4 = C4, and 
X 3 = {4/3,-10,8/3,1/3} to the case 3) C 4 = 0. We use the values of effective range 
parameters, 



a c = -7.82 fm , r = 2.78 fm , P ~ 0.022. 



(39) 



We can also have explicit expressions for the LECs C (fi), C 2 (/i), C 4 (/i) and C^fj) from 
Eqs. (34), (35) and (31) with the constraints for C 4 and C4. 

In Fig. 5 we plot our result of the strong scattering length a(/i) as a function of the 
scale parameter \i. In the left panel, we plot three curves for the strong scattering length 
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Figure 5: Strong scattering length a(/x) [fm] in functions of the scale parameter \i [MeV]. 
In the left panel, a(p) is plotted by the dotted curve, the dashed curve, and the full curve, 
respectively, for up to LO, NLO, and NNLO. In the right panel, a(/x) calculated up to 
NNLO are plotted for the three different constraints for C4 and C4, which are explained 
in the text. 



a(/i) up to LO, NLO, and NNLO with the constraint C4 = (the case 1). We find that 
the NLO correction significantly improves the estimation of a(/i), as shown by KR. 

If one looks into the details more closely, however, there is a quantitative difference 
in the results of LO and NLO between the MS and MS schemes. The value of the LO 
scattering length ai,o{p) at /i = m w in the MS scheme, which is obtained from Eq. (36), is 
a i!o = m Tr) — —30.72 fm. The LO contributions to a(/x) can be divided into three terms; 
l/a c , the term proportional to a log function and the remaining ones proportional to aM. 
Evaluating each contribution, we obtain l/a c = -0.1279, aMln (^jj = 0.0807 and 
«M(1 — Ce) = 0.0147 in units of fm -1 in the MS scheme. There is a strong cancellation 
between l/ac an d the log term which has the order of aM. Consequently l/a(/x = m n ) 
becomes a small value, making its inverse large. In the case of the MS scheme, the 
cancellation is stronger, having the log term aMln ( ^^2 ) = 0.1247 and aM(— 3C E + 
2)/2 = 0.0047 in units of far 1 . The cancellation of l/ac and the terms proportional to 
aM makes the value of l/a(/i = m n ) two orders of magnitude smaller than l/ac- As a 
result, one gets an unrealistically huge scattering length, af 5 (/i = m n ) = 738.62 fm. The 
strong dependence on the renormalization schemes of the LO contribution to a(fi) makes 
the EFT result somehow arbitrary. 

The NLO contribution, Eq. (37) can be divided into terms linear in aM and those 
proportional to (aM) 2 . The term linear in aM is comparable in magnitude with the LO 
contribution because of the cancellation in LO, as discussed above. More precisely, we 
have 1/afS = -0.0325 and -aMr fi/2 = -0.0343 in units of fm -1 . On the other hand, 
the numerical value of the contribution proportional to (aM) 2 is 0.0015 in units of fm -1 , 
which is about 5% of the terms linear in aM. 
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NLO 


NNLO-1 


NNLO-2 


NNLO-3 


a(n) 


-14.98 


-15.11 


-15.18 


-14.62 


r(/i) 




2.73 


2.78 


2.82 



Table 1: Numerical estimations (in units of fm) of scattering length a(/x) and effective 
range r(fx) up to NLO and NNLO without Coulomb effect at fx — 140 MeV. 



The NNLO contribution is very small, as can be seen from the left and right panels in 
Fig. 5 and Table 1. The reason can be easily found from the expressions for the NNLO 
terms in Eq. (38). These terms are proportional to (aM) 2 . We observed in NLO that 
the (aM) 2 t erm is smaller than the aM order term by an order of magnitude. The 
magnitude of (aM) 2 terms in NNLO ranges from about 20% to 300% of (aM) 2 terms 
in NLO, depending on the choice of the assumptions on C 4 and C 4 . Consequently, the 
NNLO correction to l/a(/i) is about 1 ~ 6% of the contributions up to NLO, depending 
on the constraints of C 4 and C 4 . 

In Table 1 we show the estimated values of the strong scattering length a(fi) and 
effective range r(fj) at \i — m n . 8 The NNLO term itself varies by an order of magnitude 
depending on the choice of the constraints on C 4 and C 4 . However, as discussed in a 
previous paragraph, its contribution to a(/x) is suppressed due to a higher order of aM 
factor. As a result, the different choice of the constraints on C 4 and C 4 affects little the 
final result, only a few percents at most. The first correction to r(fj) appears at NNLO and 
is linear in aM, whereas the NLO correction to l/a(/x) does in the aM order. Contrary 
to the case of l/a(/x) where the aM correction plays a crucial role, the aM contribution 
to r(fi) amounts to only about 2% of r . Though the aM order corrections to l/a(fj.) 
and r(n) are of the same order of magnitude, the (aM)° order contribution to l/a c is 
smaller than that of r by an order of magnitude. Consequently, we have very contrasting 
behavior of a(/x) and r(fi). 

Thus our results of the strong pp scattering length and effective range up to NNLO, 
which are estimated by employing the dimensional regularization and the MS and PDS 
schemes at /i — m n , can be summarized as 

a(/i = 77v) = -14.9 ±0.3 fm, (40) 
r(n = m n ) = 2.78 ±0.05 fm, (41) 

where the error-bars are estimated by the uncertainties due to the constraints on C 4 and 
C4, which could play a similar role to the model dependence in deducing the values of the 
strong scattering length a pp and effective range r 0iPP in the potential model calculations. 

7. Discussion and conclusions 

8 We find a minimum point for a(p) at fi ~ 2/r ~ 142 MeV, which is very close to the pion mass, 
/i = to^. 
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In this work, we calculated the S'-wave pp scattering amplitude up to NNLO in the 
framework of the pionless EFT. The loop functions were calculated by using the dimen- 
sional regularization with the MS and PDS schemes. After fixing the LECs by using the 
effective range parameters, we estimated the strong scattering length a(/x) and the strong 
effective range r(/x) as functions of fi. The LO contributions to l/a(fi) are composed of 
1/ac and the terms depending on aM arising from the loop diagrams. The smallness 
of 1/ac makes it comparable in magnitude to the aM terms in the same order. Due to 
the opposite signs of 1/ac and the aM terms, furthermore, there is a strong cancellation 
among them and thus it makes the LO result for l/a(/i) suppressed and sensitive to the 
renormalization schemes. The NLO correction, expanded in powers of aM, begins with 
the linear order of aM. The linear aM order correction to a(fi) is of the same order 
of magnitude as the aM terms in LO, and thus makes the NLO contribution crucial in 
both of the MS and MS schemes. The higher aM order terms in NLO, e.g., the terms 
proportional to (aM) 2 are suppressed to a few percents of the leading contribution, so 
they can be regarded as a perturbative corrections to both a(/x) and r(/i). The NNLO 
terms give us only a fairly minor correction to the results up to NLO. The reason is partly 
attributed to the additional order counting of the NNLO terms in powers of aM: The 
aM order corrections in NNLO begin with (aM) 2 . Similar to the (aM) 2 contribution 
in NLO, the terms in NNLO produces small corrections to the results. In conclusion, we 
can say that our investigation reveals both bright and shadowy aspects of studying the 
strong pp scattering length in EFT. Convergence from NLO to NNLO is satisfactory, but 
the LO and NLO results are significantly dependent on the renormalization schemes. 

Though the quantities of the strong scattering length and effective range from the 
pp scattering could be regarded as physical quantities, it is unlikely that they can be 
determined unambiguously without the subtraction scheme and renormalization scale 
dependence within the present framework of EFT. Similar arguments can be found in 
Refs. [30, 31]. Nevertheless, the strong pp scattering length and effective range are im- 
portant ingredients for better understanding of the isospin nature of the NN interaction. 
The problem of the strong pp scattering length may have to be approached at various 
levels, from "first principle calculations" like lattice QCD to more complex systems in 
which a(ji) (or equivalently Co(//)) plays non-trivial roles. 
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Appendix A: Amplitudes in NNLO 

In this appendix we present expressions of each of the amplitudes in NNLO in terms 
of functions, ^0,2,4 and Jo,2,22,4- Detailed calculations of the ip and J functions are given 
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in Appendix B. From the diagram (a) in Fig. 2, we have 



T. 



(4,o) 



n=0 

T^GoJoi (2^3 <^ ^ l9>y (2^ 



'(+)t>i./.(+)\ 



^ +) ^i4 +) ) 



i CqCI 



"(^0^2 + ^2 Jo) ^ 



41-CqJo 

From the diagrams (b) and (c) in Fig. 2 we have 



(4,ft,c) 



sc 



v(+)f>/v(+) 



n=0 



n=0 m=0 



m=0 



v(+)f> 



(+)\ 



" (1 - C J Q )* 7 (2^(27F (g|Gc Wc ^ 



x 



(2+J 3 



(k\G^V 2 \^ ] ) + (^>\V 2 G^\k) 



Hi 



C C 2 ?poJoJ 2 



(1 - C Jo) 2 
From the diagram (d) in Fig. 2, we have 



r (4,c0 

■sc 



<4 iE(wg- J )V Gg-v 2 Gg- ; Eiwr^v^^ EIwS'tk, 

Z=0 m=0 n=0 



(1 - C J ) 



d 3 tf d 3 q 



^|G^V 2 4 +) |g) 



(2tt) 3 (2tt) 3 



(I-C0J0) 3 ° 2 



2 7 2 



From the diagram (e) in Fig. 3 we have 



T, 



(4,e) 



-sc ~ ^<f ] \V 2 G£>V 2 \^ , { 
From the diagrams (f) and (g) in Fig. 3 we have 



*(+)Y> L/,(+)\ 



°2 / ,2 



(^22 + ^0 + 2^2^2). 



(+)> 



y(+)t> 



,(+)\ 



n=0 



+<4 _, lWc ^( W)^ |4 

?1=0 



( + )\ 



1 QA^o 

2 1 — CqJq 



(ij J 2 + 2^) 2 j j 2 + i)oJ J 22 ) . 
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From the diagram (h) in Fig. 3 we have 



',(-) 



,(+)\ 



m=0 
2/nf2„/,2 



n=0 



1 C=fC|^g 



4(1 -Co Jo) 2 
From the diagram (i) in Fig. 4 we have 



(3JqJ 2 + Jo J22) 



■sc 



(+h 



1 1 ~ 



(47) 



(48) 



From the diagrams (j) and (k) in Fig. 4 we have 



n=0 ra=0 

1 Co^o 



2 1 — Co Jo 
From the diagram (1) in Fig. 4 we have 



[2C 4 ^2 J 2 + C 4 (^ J 4 + 4>iJo) 



(49) 



n(4,0 



(+)t>/A(+) 



,(+)\ 



m=0 
2„/,2 



n=0 



1 « 



2 (1 - C Jo) 2 



C4J2 ~\~ C4J0 J 4 



(50) 



Appendix B: Loop functions 

In this appendix, we present ip functions (tp , ip 2 , ^4) an d J functions (J , J2, J22, and 
J4) employing dimensional regularization and power divergent regularization scheme [7, 8] . 
We first show the calculations of the ipo, 4>2, "04 functions in Eqs. (16), (19), (24). 

1. ipQ-. The Fourier transformation of the Coulomb wavefunction * s 



-ik-r 



with 



4 ±} ( f ) = E( 2Z + 1 )^/ ±j (H^(cos^), 



:i D (±)i 



(51) 



(52) 



where cos# = p-r. One has the relation, k-f— Ax[cos#cos0 + sm0sm0cos(</> — 0)], where 
r and are represented by (r, 6>, 0) and (/c, 0, 0), respectively. Now we choose = and 
then have 



2tt 



d(j)e -ikr sine sine cos * = 2tt J (-A;r sin sin 0) , 
16 



(53) 



where J n is a Bessel function and we have used the Bessel's first integral, J n (z) 
1 

2-iri" 



It d(pe izcos<t >e inil) . Using the relations, 



f d9 sin 6 Pi (cos 9) J ( - kr sin 9 sin 6)e~ lkr cos 9 cos 9 = i l \ ^-Pi (cos 9)J l+1 (-kr) , (54) 



Ji(—z) = (—l) l Ji(z), and ^(z) = ^Jj^J l+ i(z) where Eq. (54) is obtained from Eq. (15) in 
Ref. [32], we have 

oo poo 

i)f\k) = 4tt^(2Z + 1)P,(cos0) / dr^Rf) {prater). (55) 

Now we calculate by the dimensional regularization. The angular integration will 
pick up the / = part of the wavefunction, thus we have 

P OO P oo 

- / drr 2 R { +) (pr) / dkk d ~ 2 jo{kr) 
1 Jo Jo 



2y 
= 47T 



2/ (2vr) d 

(2*) 



= (2nf 2 (0"" | o °° drr^^ipr) f dpp^J^p) . (56) 

Using the relation J °° dtt^ 1 J v (t) = T ^~^ +v) ^ T (^), we have 

dpp^J h (p) = (57) 

Furthermore, fromEq. (6.64) of Ref. [24] we have i?^(pr) = e^C^F^l+ir], 2; —2ipr)e tpr , 
where iFi(a;b;z) is the confluent hypergeometric function (or Kummer's function of 
the first kind). Using the relation, J °° e~ t t b ~ 1 iF 1 (a, c; tz) = T(b) 2 F 1 (a, b, c; z), where 
2F1 (a, b; c; z) is the first hypergeometric function, we have 

/ e ipr r 3 - d iFi(l + ir}, 2, -2ipr) = T(4 - d)(-ip) d - 4 2 Fi(l + irj, 4 - d, 2; 2) , (58) 
Jo 



and thus 



^0 = (2vr) 3 / 2 (0"" ^ ie ^C7 l? r(4 - d)(-ip)-"Vi(l + «J, 4 - d, 2; 2)2^£p) 



.(59) 



There are no poles at d = 3 and 4 in Eq. (59). Using the relation 2-^1(1 + irji 0> 2; 2) = 1 
and fi d = 27r d / 2 /r(d/2), we have 

^0 = e mo C v . (60) 
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i-d 



2. i; 2 : 

= (2nf 2 

= (27T) 3 / 2 



For d = 4 we have 



(27T) 



(2vr) £ 



-/ drr l - d R { +\ W ) dpp d -*Ji(p) 
Jo Jo 2 



(2tt) £ 



je^C^-i^^aF^l + i V , 2 - d, 2; 2)2 d " 2 



^r(^)r(4-d) 



3-d 



r(¥) ' 

(61) 



^ 2 = e^C, ^ - -a 2 M 2 j , 
where we have used the relation 2-^1(1 + if), —2, 2, 2) = | — |r^ 2 . For d = 3 we have 



(d=3) 



-e^CnaMfi- + 



3-d 

where we have used the relation 2-^1(1 + if), — 1, 2; 2) = — irj, and thus we have 

1 



^2 = e"°C„ 



p 2 - «M|i - -(aM) 2 



(62) 



(63) 



(64) 



3. 

V>4 = 



J y (2^^ (A;)A; 



(27T 



,3/2 



/; \4-d Q r oo 



2/ (27r) d - 1 7o 



- / drr-^i^+V) / dpp d+ *Ji(p) 
1 Jo Jo 2 



(27T) 3 / 2 f £ 



ij ^£r^(-*) + -42; 2) 4(1 _, ) ^_^ ) r r y ) ) 



For d = 4 we have 



^4 = e""C„ ( P 4 - ^ 2 MV + ^a 4 M 4 ) , 



(65) 



(66) 



where we have used the relation 2-^1 (1 + if), —4, 2; 2) = -^(3 — IO77 2 + Iff). For d = 3 we 
have 



4 ... 



< =a = --e iao C v aMp [p z - -a 2 M 2 
3 



3-d 



+ 



(67) 
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where we have used the relation 2-Pi(l + if], —3, 2; 2) = \r}{— 2 + r? 2 ). Thus we have 

^4 = e^jp 4 - [^ + ^(aM) 2 ]p 2 + l(aM)V+^(aM) 4 }. (68) 



Now we calculate loop functions J , J2, ^22 and J 4 in Eqs. (17), (20), (22), (25) by 
using the results of the ip functions obtained above. 



4. J : 

The function Jq(p) is given by [7] 



d z l 2ttt7(0 



(2tt) 3 e 2 -"(0 - 1 p 2 - Z 2 + ie ' 
where / = |/|. We now separate J into two parts as [7] 

Up) = Jt + 4 in , 



where 



(69) 

(70) 
(71) 

(2tt) 3 eWO - 1 Z 2 p 2 - / 2 + ie ' ^ 

As J(( m is already calculated in Ref. [7], by changing the parameter x = 2txt](1) and 
using the relation 



jdiv 
J 



jfin 



-uj 



d 3 l 2tt?7(Z) 1 



M 



r or, 
J (2tt 



(2tt) 3 e 2 ^W - 1 Z 2 ' 
d 3 f 2ttt7(0 1 V 2 



f 

Jo 



dx 



x 



(e x - l)(x 2 + a 2 ) 2 



In 



1 \ 7T 

2tt/ ~ a 



2tt 



where ^ is the logarithmic derivative of the T-function, we have 



J< 



fin 



aM 2 



aM 2 



M 



Hv) - c;-(i P ) , 



All 47T V 4:7l 

where 77 = aM/(2p), H(r)) = ij}{irj) + ^- — ln(ir]), and h(rj) = ReH(rj). 
Next we calculate the divergence part Jq iv in d = 4 — 2e dimension 

4-d , d d-iy> 2it<q(q) 1 



jdiv 

J 



-M 



(27r) d - 1 e 2 ^(i) - 1 g 2 • 
Changing the variable x = 2nr](q) = naM/q, we have 



div 



/..\4-d n <d-l)/'z 

\2J (OttV-it J Jo 



= —M 



(2n)^T (i 2 
//A 4-d 2 7r( (i - 1 )/ 2 



3-d 



e x - 1 



(73) 



(74) 



(75) 



TY (a7rM)^ 3 r(4-rf)C(4-rf), 



2/ (2 7 r) d - 1 r(^i 
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(76) 



where we have used the relation Q d = 2n d / 2 /T(d/2) and ((z) is the Riemann's zeta 
function. For d = 4 — 2e we have 



j; 



div 



aM 2 

~87T 



7- 3 ^ + 2 + ln fe 



(77) 



We also consider the pole for d — 3, known as the power divergence subtraction (PDS) 
scheme pole. Using the relation lim^x £(s) — = 7, we have the pole at 3-dimension 



J ° 47T3-rf + 

and thus we include the PDS counter term and have 



dlv M aM 2 



(78) 



(79) 



5. Jo 



o? 3 g ^2(5)^0(5) 



(2tt) 3 (2tt) 3 
Using the result of tp2 in Eq. (64), we get 



■Up) = 



p 2 - (uxM - - (aM) 2 



(2tt) 3 p 2 — a 2 + ie 



J (p) - AJ 2 , 



where 



AJ 2 = M 



4-d f - rid— I 



d d - x k 2nr](k) 



.2 J J (2tt) o! - 1 e 2m '( fc ) - 1 
For d = 4 we have 



M 



2) 4 "9^ (7raM)<i " ir(2 - rf)c(2 - d)m 



(80) 



(81) 



AJ 2 = -7ra 3 M 4 C'(-2) , 
where C'( -2 ) = -0.0304 • • •. For d = 3 we obtain 



AJ. 



(d=3) 



— %a 2 M 3 u + 

48 P 3-d 



and by including the PDS counter term we have 



AJ 2 = ^(aM)V+^f(«M) 3 C'(-2). 



(83) 



(84) 



(85) 



6. J, 



22 



J22 — 



d 3 q dV j 2/ j lM+) ,_, _e _ % , f d 3 q Ml)^) 



(27r) 3 (2vr) 3 ^ Wl ° 7 (2tt) 3 p 2 - g 2 + te 

(p 4 - 2Ap 2 + A 2 ) J - {p 2 - 2A)AJ 2 - AJ 22 , 



(86) 
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where A = jj,aM + \{aM) 2 , and 



AJ 22 = M(| 



(2tt) c 



= M, |) 4 "(|p L T^ M )' i+lr (- rf )^- rf )- 



For d = 4 we have 



(87) 



AJ 22 = -fVMT(-4), 



where C'(~4) = 0.00798 • • •. For d = 3 we get 



AJ^= 3) 



'22 



and thus we obtain 



AJ, 



22 



.-L^M 5 //— ^ + 
2880 P 3 - d 



■S( a M)V + ?(«M)V(-4). 
2880 v ; P 48 v ; v ; 



(88) 



(89) 



(90) 



7. Ja 



d 3 g Mq)^o(q) 



(2tt) 3 (27 

Using the relation for i\) A in Eq. (68), we have 



(2tt) 3 - q 2 + ie 



J 4 



P 4 - 



4 5 
-al/i + -(aM) 2 
3 6 

p 2 -^aM/i-^(aM) 2 



p 2 + I(aM) 3 /x + l(aM) 4 } J 



A J 2 - A J 
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(91) 



(92) 
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